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ABSTRACT 

The    spectrum  of   the   mono-energetic   neutron  transport 
operator  for  a  homogeneous    sphere   where    the   neutrons    are 
scattered  with  spherical   symmetry  in   the    laboratory 
system  is    considered.      The   operator  may  be  written: 

Af  ^  .  <*fj^i?'      .   c-  f  (x,y)    +  £  0    (r) 

r 
0(r)    =    f       f(z,Vr2-z2)    dz 


7-r 

2    2    2    2 
in  the  region   r  =x  +y  <  a   ,  0  <  y  with  the  boundary 

r~z  2~ 

condition  f(-Va  -y  ,  y)  =0  and  with   cr  and   c   real, 
positive  constants. 

The  characteristic  value  problem  Af  =  Xf   is  to 
be  analyzed.   An  associated  integral  transformation 

a  | r+s | 

V-  [  0{s)  ds  [        e'Xw^r 

4>         J|r-s| 

arises  in  a  natural  way.   For  real  X   the  operator  L. 
is  hermitian.   The  characteristic  value  problem 

will,  with  the  relation 
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f(x,y)  =  %   f       ft(\fzZ+jZ)    e" 


-X (x-z) 

dz 


'a  -y 

yield  solutions  of  the  problem  when   p  =  2/c  . 

The  real  spectrum  of  A   is  analyzed  by  investigation 
of  the  integral  transformation  K   for  real  X  „   In 
this  case  it  is  shown  that  its  characteristic  values 
p. (X)   i  =  1,2,...   are  positive  and  depend  continuously 
on  X  o   Comparison  transformations  show  that 

lim   p.(X)  =  0  and    lim    P*(X)  =  oo  . 
X  — >  oo  1  X  — >  -  oo 

This  proves  the  existence  of  an  infinite  number  of  real 

characteristic  values  of  A  .   These  characteristic  values 

are  shown  to  have  index  one. 

The  proof  is  based  on  an  extension  of  methods  used 

by  Lehner  and  Wing  [jjo 
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1.   Introduction 

The  neutron  density  at  time   t   in  a  body  may  be 
described  by  a  non-negative  function  \|/(x,v,  t)   where 
x  =  (x,,x2,x  )   and   v  =  (v1, v2,v^)  .   This  function 
satisfies  the  linear  Boltzmann  equation 

dl  =   -v    .    gradx  \|f    -    cr(x,v)    \|>   +   J     K(x,  v,  v'  )\|/  (x,  v' ,  t)dv«       (1.1) 

for  x   in  a  domain  X  ,  v   in  a  domain   V   and  positive 

t  .   Usually  X   is  assumed  to  be  convex;  in  this  case  a 

standard  boundary  condition  is   \|/(x,v,  t)  =  0  for  x  on 

the  boundary  of  X   and   v.n  <  0  ,  where   n   is  the 

outward  normal  to  the  domain  X   (i.e.  no  neutrons  enter 

X   from  outside   X) .   The  domain   V   is  bounded.   The 

non-negative  function  cr(x,  v)   is  the  macroscopic  removal 

cross-section  multiplied  by  v   for  neutrons  at  x  with 

speed  v  .   The  non-negative  kernel  K(x, v, v' )   is  a 

function  or  distribution  and  the  integral  transformation 

describes  the  production  of  neutrons  by  fission  or 

scattering  at   x  with  velocity  v  „ 

The  initial  value  problem  is  to  find   \|/(x,  v,  t) 

for   t  >  0  which  assumes  an  initial  density   ty(x, v, 0) 
as   t  — >  0+  . 


-  5  - 


K.  Jorgens  has  shown  under  rather  general  conditions 
that  this  is  a  well-posed  initial  value  problem  for  X  x  V 
boundede   The  existence  of  a  strongly  continuous  semi-group 
operator,   E(t)  ,  which  solves  the  initial  value  problem 
is  demonstrated.   In  some  cases,  additional  assumptions 
assure  that   E(t)   is  compact  after  some  finite  time  which 
depends  on  the  domain  X  x  V  9    the  character  of  the 
scattering  operator,,  and  the  underlying  Banach  space  on 
which  the  problem  is  posed0   In  such  a  case  the  spectrum 
of  A   is  discrete,  confined  to  a  left  half  plane  and 
finite  in  every  finite  vertical  strip  ([1]  Th„  5°l)e 
Furthermore,,  an  asymptotic  expansion  of  E{t),   the  semi- 
group operator,  in  terms  of  the  projection  operators 
associated  with  the  generalized  characteristic  functions 
of  A   exists   CEU  Th,  5.2). 

In  an  earlier  series  of  papers  Lehner  and  Wing 

([2],  [3],  [I4J )  have  completely  analyzed  the  spectrum  of 

the  transport  operator  for  a  simple  case,,  the  infinite 

slab  of  finite  width  where  the  neutrons  have  only  one 

speed  and  are  scattered  with  spherical  symmetry  in  the 

laboratory  system.   This  problem  may  be  written 

1 
|4|  Cz,ii,t)  =  -  p.  j-f  -  er*  +  f  C   *(z,(i«st)  dji'   (1.2) 


\|<(a,u.,  t)    =   0  for     [i  <   0 

ty(-a,  y,,  t)    =   0  for     u-  >   0 

\|/  (  z,  |x,  0) 


is   given   on 


X:       [  z  |    <   a   5     |u.|    <   1 

V.  :      v   ,    a   positive    constant 

\|f    6   L2(  |z|    <  a   ;    \\i\    <   1) 

Here   0s  and   c   are  positive  constants.   The  domain  X  xV 
is  unbounded  and  Jorgens'  conclusions  do  not  apply.   The 
point  spectrum  of  the  operator  is  shown  to  be  finite  but 
not  empty,  the  continuous  spectrum  fills  a  left  half 
plane,  and  the  residual  spectrum  is  empty.   Lehner  also 
partially  analyzed  the  same  problem  where   X   is  a  finite 
sphere  ([ij.]).   Here  he  determines  in  accordance  with 
Jorgens  that  the  operator  has  a  discrete  spectrum  with 
-OD  the  only  possible  accumulation  point. 
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2.   The  Mono-energetic  Neutron  Transport  Operator  for 

the  Homogeneous  Sphere,, 

With  the  following  simplifying  specializations 
(i)   X   is  the  domain   r  =  |jcj   <  a   , 
(ii)   V   is  the  domain   jyj  =  const,  >  0  , 
(iii)  cr   (xsv)   is  a  positive  constant, 
(iv)   Scattering  is  spherically  symmetric  in  the 
laboratory  system, 

the  problem  (1»1)  becomes 

+  |     \|f(r,ix',t)  d»i' 

with  boundary  condition 

\l»(a,  (j.,  t)  =  0        for  [i   <   0 
and  the  initial  condition 

\|/(r,  m  0)   prescribed  on  X  x  V  . 

Here   R   is  the  distance  from  the  center  of  the  sphere, 
\x      is  the  cosine  of  the  angle  the  neutron  beam  makes 
with  the  radius  vector,  and   o^  and   c   are  real,  positive 
constants  s 

It  will  be  convenient  to  perform  a  change  of 
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variables  before  rendering  the  problem  more  specific. 
Let: 

x  =  r  [i. 


y  =  rV  1-p. 

\|r(r,n,t)  =  f(x,y,t)   ; 


then 


^f(^ty>t}  =  Af  (2.2) 


^'U'  (Tf   +  2%  i"   r(xls\/7^f",t)  dxx 


with  the  initial  data   f(x,  y, 0)   prescribed.   The  operator 
A   is  defined  over  a  Hilbert  space   H   of  square  integrable 

functions  with  the  inner  product 

s/"~2   2 
a       J  &   -j  


(f,g)  =  kn    \  y  dy  \   f(x,y)  g(x,y)  dx  . 

At  first  glance  this  inner  product  may  seem  unusual,  but 
it  is  merely  the  transformation  of  the  natural  inner  product 
for  the  operator  in  terms  of  the  original  independent 
variables . 

-+1 


(\|>1,*2)  =  lpt  V  r2  dr  I    du-  ^ 


The  domain,   D(A)  ,  of  the  operator  in  (2.2)  consists  of 
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those  functions   f  which  satisfys 

(i)   f(x,y)   is  defined  for  r2=x2+y2  <  a2   and  y  >  0  , 
(ii)   f  e  H  , 
(iii)   Af  e  H  , 


In  particular,   Af  is  absolutely  continuous  in  x  for 

all  y  . 

This  is  a  well-posed  initial  value  problem;  the 
operator  has  a  discrete  spectrum  with  ~oo  as  the  only 
possible  accumulation  point  ([l]  and  [1+])  . 

If  f(x,y)   is  in  D(A)   and 


(*  +  o)    f(x,y)  =  -&  +  &[        f( 

U-r 


z,\/r2-z2)  dz 


then  eat//v  f(x,y)   satisfies  (2.2)  with  initial  data 
given  to  be   f(x,  y)   and   f   is  a  characteristic  function 
of  A  with  characteristic  value   a/v  , 
Abbreviate  X  =  a/v  +  C    and  define 

rF 


0(r)    =  I    f(z,  v/r2=z2)dz  .  (2.3) 


-r 


(Notice  that  0     is  an  odd  function.)   Rewrite 
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Since      f   e   D(A) 

^(\/x2+y2)    -S ^ 


^7  ^W 


and   from    (2.3) 


J2f(r)    =    f     f(x2,  \JtZ-x\) 


dx~ 


I  f  °"X*2  *»2f2,,   ,  ^<^!^I> ,—^7  tai 


J-Va   -r  +x|  Vxi+ 


2 


The  change  of  dummy  variables 

2    2    2    2 
s   =  x^  +  r  -  x| 

w  =  Xp  -  x1        yields 

0(r)  =  |  1^0       where  (2.I4.) 

L^  =  [    0(a)    ds         e-Xw^ 
J0  J  I  r-s  I 


This  integral  equation  will  be  the  focus  of  attention. 
Theorem  2.1.   The  necessary  and  sufficient  condition  for 
f(x,y)   to  be  a  characteristic  function  of  the  transport 
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equation  is  that  its  associated  0  be  a  characteristic 
function  of  the  integral  equation   p0  -   L.  0  with 
characteristic  value   p  =  2/c  „ 

Proof0   The  derivation  of  the  integral  equation  has  proved 
sufficiency,,   The  necessity  follows  by  reversing  the 
computation,,   That   0(r)   as  defined  in  (2.3)   is  in 
LpCO^a)   is  a  consequence  of  f  e  D(A)   and 


(0,0)  =  I   dr 
^0 


\/2   2. 
.  V  r  -z 


f(z,Vr  -z  )dz 
■r 


<  \   dr  2r  \ 


f(zs  Vr2-z2) 


dz 


a 


/ a2-y2 


=  2  \  y  dy 

J0      V,1  .v/a^Ty^ 


\f(xtj)\d   dx 


2k 


Cfa  f)  <  oo 


In  what  follows  it  may  be  assumed  that  0  has  norm 
one9  since  no  explicit  use  will  be  made  of  the  norm  of  the 
characteristic  functions  of  the  transport  operator  except 
that  it  is  not  zero, 

The  knowledge  of  solutions   0(r^A)   of  (2el|)  will 
provide  information  about  the  spectrum  of  the  transport 
operator,,   The  integral  equation  p0  =  L.  0  resembles  the 
usual  characteristic  value  problem,,   However,   p   is  a 
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positive  constant  dependent  on  the  nuclear  properties  of 
the  materials  in  the  sphere.   Therefore,  only  characteristic 
functions  0(r^X)      with  positive  characteristic  values,   p  , 
are  to  be  considered  and  only  those  with  characteristic 
value   2/c  will  be  the  characteristic  functions  sought,, 
A  list  of  properties  of  L,   follows. 

Theorem  2,2.   The  integral  transformation,   L   ,  defined 
for  all  complex  X   with  domain  D(L  )   taken  to  be  the 
Hilbert  space   Lp(0,a)   with  the  ordinary  inner  product, 
has  the  following  properties , 

(i)   Its  kernel  is  symmetric  but  not  hermitian 
(unless   X   is  real) , 

(ii)   It  is  compact, 

(iii)   It  is  non-degenerate  for  real   X  . 

Proof.   Property   (i)   is  obviousc   Notice  that  for  X 

real  the  kernel  is  a  positive  function,,   Therefore  the  trans< 

formation  maps  non-trivial  non-negative  functions  into 

positive  functions.   Coupled  with  the  compactness  of  the 

operator  this  will  imply  that  the  dominant  characteristic 

value  of  L   is  simple  for  real  X   ([7]). 

The  kernel  of  L    (2<>k)   may  be  written  as   |  r-s  | 

A. 

B(r, s)   where   e   is  an  arbitrarily  small  positive  number 
and   B(r, s)   is  a  bounded,  continuous  function  of  r  and 
s  ,   Therefore,  the  transformation  is  compact  ([6],  p„326). 
If  U   were  degenerate  for  real  X  ,  its  kernel  could 
be  written  as  a  finite  sum   <r —  p„0,(r)  j^TTsT  where  the 
0.      would  be  the  normalized  characteristic  functions  of 
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the  operator.   Therefore  the  trace  of  the  kernel  would 
exist  and  be  finite 0   However,  this  is  not  the  case. 

In  short   L^   isP  for  real  X  9  an  integral  trans- 
formation of  Hilbert-Schmidt  type 0 

If  0     is  extended  as  odd  to  the  interval   |r|  < 
then  for  (3  5  Re(X)  >  0  ,  (2ol+)  may  be  rewritten 

"  .)    =    f  0(s)    ds    f  e^w  ^  -    (  0(s)    ds    f  e"Xw  dw 


a 


p  jflir;    =    |      Jfl^sj    as    1  e    -      —  -  p\s)    aa    1  «  — 

Oo  ^|r-s|  ^0  ^|r+s| 

fa 
=    \      0(s)    E(X|r-s|)    ds    -  0(s)    E(\|r~s|)    ds 

a 
p  0(r)    =    f       0(a)    E(XJr-sj)    ds  Re(X)    >  0  (2.5) 

which  is  the  integral  equation  associated  with  problem 
(1.2)  for  the  infinite  slab  of  width  2a   ([2]). 

Lehner  and  Wing  ([2],  [3])  showed  that  (2.5)  has 
solutions  with  positive  real   p   only  if   ^=  Im(X)  =  0 
and  p  >  0  o   In  the  case  0-   0  and   (3  >  0 

J2f(s)  E(plr-s|)  ds 

is  a  positive  definite  operator  with  positive  character- 
istic values   p,  >  Pp  >   000   and  for  each  i  ,   p.   is  a 
continuous  decreasing  function  of  (3  which  tends  to 
zero  as   (3   tends  to  infinity.   Also,   p,   approaches 


Hi- 


infinity  as   (3   approaches   0   from  above  ([2]),   This 
is  summarized  in  Figure  10   These  characteristic  functions 
may  be  taken  to  be  odd  or  even.   In  particular,  the 
dominant  eigenvalue  is  simple  and  corresponds  to  a 
characteristic  function  of  one  sign,   hence  even  ([7JS 
Th»  6o3)<»   Deleting  the  even  functions  in  order  to  obtain 
solutions  for  the  sphere  might  lead  to  a  graph  like 
Figure  2.   There  will  be  at  most  a  finite  set  of  values  of 
P  ,  and  possibly  none,  which  yield  characteristic  value 
2/c  o   The  corresponding  characteristic  functions  and 
values  of   p   for  this  integral  equation  problem  will 
yield  characteristic  functions  and  characteristic  values 
for  the  transport  equation,, 
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a> 


3 

on 
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3<>   The  Real  Spectrum  of  the  Operator 

The  integro-differential  operator  is: 


0(r)   =     f(z,  v^r2=z2)  dz 


(3.1) 


on 


X  s  r2  =  x2  +  y2  <  a2 


y  >  0  . 


The  domain  D(A)   consists  of  functions   f  in  a  Hilbert 
space   H  of  functions  square  integrable  on  X  such  that 


f  e  H  2   Af  e  H 


f(-Va2=y2py)  =  0 


where    the   inner  product   is   defined   to  be 


(f,g)    =  U.it    \     y  dy 


tt? 


^-r 


fCx,y)    g(xf y)    dx    . 


The  adjoint  operator  acts  on  the  same  Hilbert  space  and  is 


*f  =  +  m^Lll     -   (Tf(x,y)    +  &  0(r)    .       (3.2) 


^Tx 


2r 


,  * 


Its   domain     D(A   )      consists    of   functions      f      such  that 


2      2 


f  e  H  |     A"f  e  H  9      f(+\/aVsy)   -  0 


Since   o~  (A)  ,  the  spectrum  of  A  ,  is  discrete. 
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<T(A)  =  (T  (A*)   and  the  spectrum  of  A*   is  discrete 

([8],  p.  11+3)  o   Notice  further  that  if  X  e  cr  (A)      with 

characteristic  function  f(x, y)  ,  there  naturally 

corresponds  to  it  a  X  e  ^(A)   with  characteristic 

function  f U,y)  .   Thus,  the  spectrum  of  A   is  symmetric 

with  respect  to  the  real  axis  of  the  complex  plane. 

If  f(x, y)   is  a  characteristic  function  with 

characteristic  value  X  ,  then 

*a         J  r+s  | 
§  0(r)    =  0<s)  ds  1      e'Xw  ^  ,         (3.3) 

Jo        ^|r-s| 

and 

rXx-j^ 
t j*0/xf+y2)  g dx   . 

Define   Re(X)  =  P  =  -Y  °   If  only  real  values  of  X   are 

considered,  then  the  integral  operator  has  a  real 

symmetric  kernel  and  hence  only  real  characteristic 

values. 

Theorem  3.1.   The  integral  transformation  L   is  positive 

definite  for  real   (3  s 

Proof.   The  kernel  of  the  transformation  may  be  taken  to  be 

k(r,s^p)  =  k(r-s^p)  -  k(r+s|p) 


where 
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.2a 


-Sw   dw 

e  w- 


for      |q|    <   2a 

k(q;P)  =   \  VJIql 

(       0  otherwise 

If      p  >  0   ,    consider  the   fourier   cosine    transformation 
([9],    p. 314-6) 


<v>  ->/! 


2a 
cos    vz  dz    1      e 


-pw  dw 


w 


So 


and 


rg   0       e-Pw     3in  vw 
'  V  rc    \  w  v 


dw  >   0   . 


1 


,00 


k(r,s;P)    SJ=-   \         (cos   v|r-s|    -   cos   v    |  r+s  |  )    f    (v)    dv 


n  r 

~   2./—    \        sin  vr  sin  vs    f 

V*  Jo 


(v)    dv 


(0,Lp0)    =    I      0(r)    dr    (      ^TiT  ds^   (      sin   vr  sin 


vs    f    (v)    dv 
c 


Only  functions      0     which  are   odd   and   zero   outside    the 
interval      (-a, a)      need  be  considered.      Thus, 


(0.1*0) 


=  [J  f 


p"   "V  3? 

^0 


(v) 


0(s)    sin  vs   ds 


Jo 


dv  >  0   . 


If      p   <   0   ,    then     y      is   positive   and 
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De 


r+s| 

k(r,s;-Y)  =  1       e*w^  (3M 

r-s  | 

"  r+s  | 

(cosh  yw-1)  ~   + 
r-s  | 

(3.5) 

' r+s |  A  r+s | 

+  2  ('     (1-e-Y")  &+  I     e-Yw  ^L 
r-s  I  0|  r-s  I 


note  these  three  terms  by  k,  (r,s|-Y)  *  k_(r,s;-Y) 


and  k.(r,  s|-y)  respectively. 

The  first  part  of  the  proof  shows  that   k   is  the 
kernel  of  a  positive  definite  transformation.   The  same 
technique  suffices  for  k_  ,  which  is  then  also  the  kernel 
of  a  positive  definite  transformation. 

Finally 

n|r+s| 
k^r.si-Y)    -  2    \  (cosh  Yw   -   1)   ~ 

Jl r-s | 

' r+s |  v 

dw  sinh  wz   dz 

II  r-s  |         JO 


=  2    \      (cosh  z | r+s |    =   cosh  z|r-s|)   — 


i|    \     sinh  rz    sinh  zs  — 
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The  re  f o  re , 
.a 


a  a  y   a 

(  0(r)    dr  \  ^TTT  ^(r^s^-r)  ds  =  1+  f  {    0(a) 
J0         J0  ^0  T> 


sinh  zs  ds 


dz 


The  integral  transformation  L   is  a  sum  of  positive 
definite  transformations  and  hence  is  positive  definite 
for  all  real   (3  0   Its  characteristic  values   p.  (P)>PP(P)> 
satisfy  pi(P)  >  0  for  all   i   and  all  real   j3  «,   The 
characteristic  value   p,   is  simple  and  has  a  positive 
characteristic  function  ([?]*  ThQ  6„3)8 
Theorem  3o2.>   For  each  i  ,      p  (p)   Is  a  continuous 
decreasing  function  of   p  0 

Proof.   The  addition  of  a  bounded  positive  definite 
transformation  to  a  bounded  positive  definite  transforma- 
tion increases  every  characteristic  value  ([5Ji>  p0239)0 
This  induces  a  natural  partial  ordering  among  bounded 
positive  definite  operators 0   In  the  sequel   A,  >  A 
will  mean  that   A-,p  A^   and  A-i-Ap   are  bounded^  positive 
semi-definite  operators. 

The  proof  is  simplified  by  considering  two  separate 
cases  o 

(i)  o  <  p1  <  p2 


Consider 


r+s 


k(r,  sjp1)  -  k(r,  s.^p2)  = 


(e 


•  P-jV 


r-s  | 


'    w 


and  apply  the  fourier  cosine  transform  to  obtain 
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8  C 
k(rss^P1)  -  k(rss;p?)  =  ~  I   sin  vr  sin  vs  d,(v)  dv 


Jb 


,2a 


-p  w   -p  w  gin  yw 
1    =  \         °e     ' 


vw 


dw 


'0 


-p,w  -p  w 
Since   (e    -e     )/w   is  monotone  decreasing  for 

0  <  p,  <  P2   and  w  >  0  ,  it  follows  that   d-^v)  >  0 

Thus, 

V  f5(r)    dr  P  FTsT  (k(rss^p1)  -  k(rss;p2))  ds 
0O       Oo 


dl(v) 
J0 


J#(s)  sin  vs  dx 


;0 


dv  >  0 


fii)  -r1   <   -Y2  1  0 
From  (3.1;) 


rYl 

:(r,s|-Y1)-k(r,s;-Y2)  =  ^  V    si 


dz 

inh  rz  sinh  zs  —  + 

z 


Y2 
r+s| 

r-S  | 


(e-r2»  .  e-rx»  *, 

W 


The  first  term  on  the  right  hand  side  is  clearly  positive 
definite.   The  second  term  is  exactly  like  the  kernel  of 
part  (i)  and  so  this  part  is  proved. 
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That  each   p,(P)   is  a  continuous  function  of   p 
is  a  consequence  of  the   L   continuity  of  k(r,  s^(3) 

in   p   ([51,  p.239). 

The  next  step  is  to  show  that   pi(p)   tends  to   oo 
as   p   tends  to   -  oo  ,   First  observe  that 


k,  (rss£-Y)  =  k-    1   sinh  rz  sinh  zs  — 

^0 


=  k    \      sinh  %■  rw   sinh  ^  ws 

a  aw 


dw 
smn  j-  rw   sinn  ■*■  ws 
a  a 

J0 
Define 


[T1(Y)]2jZT=    \     0(s)    ^(r.s^-r)    ds  Y  >  0  (3-6) 

where      T    (y)      is    the  unique   positive   definite   square    root 
of   the   right  member  of    (3.6)    (iSh    Po265),      Define 

[T0(Y)]2#=     \     0(s)    k0(rr3|-Y)    ds  (3.7) 

^0 


v  (r0ssY)  =  -  I   sinh  ^  rw  sinh  ^  ws  dw  » 
0    "  *   —  a  \       a         a 


Lemma  3<>3»   In  the  natural  partial  ordering  of  bounded 
positive  definite  operators   [Tq(y)]2  <  [T-^y)]2   . 
Proof,   The  difference  of  the  squares  of  the  operators 
is  positive,,  since 


2k 


l2     rn2 


Pa 


«a 


(0,(Tj-Tg)0)    =  1M    0( 

Jo 


r)    dr   V   j#(s)    ds    \    (—  -— )    sinh  %■  rw   sinh  £  ws    dw 


Jo 


^0 


<H> 


•sa 


0(s)    sinh  ^  ws    ds 


'0 


dw  >  0 


for     y  >  0 


By  inspection,      ^M^      ma^  ^e  written 


TnCy)    0  =  —    \     0(a)    sinh  *•  rs    ds  for  y  >  0    „      (3.8) 


'0 


Lemma  3<,i|,   The  transformation  ^(y)   has  the  following 
properties; 

(i)   It  is  the  unique,  bounded,  symmetric,  positive 

o 

definite  square  root  of   [T  (y)]   „ 

(ii)   It  is  a  compact  transformation  of  Lp(0,a)   into 
itself, 

(iii)   It  is  a  non-degenerate  transformation. 
Proofo   The  transformation  is  clearly  bounded  and 
symmetric   The  kernel  has  a  uniformly  convergent  Taylor's 
expansion 


2       v       2 

—  sinh  -^   rs  =  — 
a     ^ 


which  is  of  the  form 


k  (r)  k  (s 
n     n 


x2n+l 
rs  ; 


Therefore, 


(3.9) 


the  transformation  is  positive.   Direct  computation  shows 

2 

that    (3„8)    is    a    square    root   of      T        and  uniqueness    is 

immediate    ([5]j>    p. 265)  . 


25 


Compactness  follows  directly  from  the  boundedness 
of  the  kernel  ([6], p.  326). 

If  T   were  of  finite  rank,  then  there  would  be  a 
non-zero  odd  polynomial  p(s)   such  that 


p(s)  sinh  -L   rs  ds  •  ■  0 

9.  ^~ 


However, 


ds  = 


=  „  *I1    f   P(t-a)  e  a    dt 

This  last  expression  is  a  Laplace  transform.   By  the 
uniqueness  theorem  for  the  Laplace  transformation  ([12], 
p.  35)  the  expression  can  vanish  only  if  p  =  0  almost 
everywhere  in  the  interval  of  integration.   Therefore 
T0  is  not  degenerate. 

All  the  characteristic  values  of   T~   have  index  one, 
because   T0   is  symmetric. 

Attention  is  now  directed  at  the  behavior  of  the 
characteristic  values   e. (y)   of   T 0(y)   for  y  >  0  „ 
Lemma  3.5.   For  each  i  =  1,2,...   the  characteristic 
value   e.Cy)   is  a  positive,  monotone  increasing  function 
of  y  which  tends  to   oo  as   y   tends  to   oo  . 
Proof.   That   e. (y)  >  0   is  shown  in  the  preceding  lemma0 
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Consider  T  (ky)  -  kTQ(y)  .   This  difference  has 


the  kernel 


—  (sinh  ^  rs  -  k  sinh  £  rs)  (3.10) 

i  a  a. 

s/a 


oo  „_2n         2n+1 


,— -  >   (2n+l)i  K    a  ; 
v/a    0 


As  for  the  original  kernel,  this  function  has  a  uniformly 
convergent  Taylor's  series  in  powers  of   rs   with 
positive  coefficients.   Since  this  is  also  the  kernel  of 
a  positive  definite  integral  transformation,  for  k  >  1 


e^ky)  -  ke^y)  >  0  (3.H) 


Then  for  y  >  Y0  >  °   it  follows  that 


e, (Y)  >  e.(Yn)  *-  >   0  (3.12) 

1       i   u  Yq 


Therefore   e. (y)  tends  to  infinity  as   y   tends  to 

infinity. 

Theorem  3.6.   The  characteristic  values   p.  ((3)   of   L„ 

tend  to   +  oo  as  [3   tends  to   -  oo  .   Indeed,  for 


P  <  P0  <  0 


s 


p.(p)  >  Ee1(r0)  f-]2   . 


Proof.   By  the  preceding  lemma,  for  some  fixed  positive 
y  =  -  p   ,  the  characteristic  values  of   Tr/Y^   satisfy 
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ro 

<D 

w 

3 
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e±(r)   >   ei(r0)  if-  for  Y  >  Y0  . 

The  characteristic  values  of   [Tq(y)]    are   [e. (y)]   . 

Since   Lp  >  [T-^y)]2  >  [T0(y)]2   for  p  <  0  ,  the  desired 
inequality  is  valid  and  the  theorem  is  proved. 

The  behavior  of  the  characteristic  values  of   p.  (p) 
is  summarized  in  Figure  3»   The  graph  may  be  somewhat 
idealized,   For  example,  it  may  be  possible  that  some  curves 
intersect  or  are  not  always  convex. 

A  brief  recapitulation  of  the  properties  of  L„ 
and  their  implications  for  the  transport  operator  follows. 

(i)   For  every  fixed   p   the  operator  has  a  countable 
infinity  of  bounded,  positive  characteristic  values  which 
accumulate  only  at  zero. 

(ii)   Every   p±((3)   is  of  finite  multiplicity,  and 
of  index  one.   In  particular   p-.(P)   is  simple  and 
corresponds  to  a  function  which  may  be  taken  to  be  positive 
except  at  the  origin. 

(iii)   Every   p. (p)   is  a  monotone  decreasing, 
continuous  function  of   p   and 


lim   p  =  0    and    lim   p  =  od 
p_>oo  x  p_i>_oo  1 

Therefore  for  every  p^(P)   there  is  a  unique   p.   such 

that   p1(pi)  =  §  >  0  0 

By  using  the  second  relation  of  (3.3)  a  characteristic 

function  of  A  ,  the  transport  operator,  may  be  found  whose 

characteristic  value  is   p.  -  o~*  . 
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(iv)      Therefore  the   transport   operator  has    a   countable 
infinity  of   real   characteristic   values   accumulating   at 
minus    infinity.      Each   of  these   characteristic   values 
corresponds   to   a   finite   number  of   characteristic   functions. 

Finally,    it   will   be    shown  that    the   real   characteristic 
values      p.      of  the    transport   operator  have   index  one. 
Lemma   3.7.      If     Af  =   pf     and      f'"'(x,  y)    =   f  (=x7y7,      then 
ff,f*)    ^  0    . 
Proof.      A   lengthy  computation  is   necessary. 


>/a2-y2 


p-  (f>f":">  =  \y d? 


c2 


/    2       2 


f(x,y)    f(-x,y)      dx 


>a 


2      2 


,*>/a  -y 


y  dy 


dx 


dp 


;0 


s/a2-y2^\/a2-y2^x 


Va2-y2 
dq 


■Mq-p)   0(  vV+y2)   0(\Iq2+y2) 


s/7^? 


vV+y2 


c2 

T 


a  w 

\      dw 
^0 


,w 


dx    \       dp    \    dq 
w        0-w        ^x 


w   e-X(q~p)    0(  v/p2-«-a2-w2)0(V/q2-Ha2-w2) 

v/q2+a2-w2 


/    2      2      2 

y/p   +a   -w 


c2 


,a 


dw 


'0 


-W 


dq    \       dp    |p-q|    w   e 

'-W         J-w 


-*lp-q|    0(j/p2+a2-w2)0(  \/q2+a2-w2 

f'  2     2     < 
p  +a  -w 


/r2 


^2 


q  +a   =w 


'0 


dw    1    dq 
'0      Jo 

.a         ^a 
dr   \    ds 


dp    (|p-q|    e-A'P-^l    +    |p+q|    e^lp+ql) 


w 


0 
a 

4 


2      2 
a   -r 


dw  0(r)    0(s)    w 


1  1     !    4k(r)-5(s)| 

fTrT "  rriir  + 


+   llTiT  +  fTTT 


1      ^k(r)i(s)| 


) 
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where 

£(r)    =     V?-aV 


then 


£    (f,f*)    =  |f   f  dr    f  ds   0(r)    JZl(s)    (e-Xlr-sl-e"Xlr+sl)       (3.13) 

Jo        Jr 

=  °£    fdr    Us    0(r)    0(a)     (e-Xlr-sl.e-X'r+sl) 

Jo      Jo 

In  the    case   where     X   =   p      is    real,    then   the    characteristic 
function     0     may  be    taken   to  be    real. 

£   Cf.f*)    =^    fdr    Tds   |Zf(r)    0(s)     (e-Pl-5!    -   .-Pl**«l) 

Jo      Jo 

If      p   >  0   j,    apply  the    fourier   cosine    transformation 

to    the   kernel   of   the   quadratic    form    (3<>13)»      This   kernel 

is    to   be  modified   exactly  as    in  Theorem  3.1. 

J  r+s  | 
1    (e-P|r-s|    _   e-P|r+s|}    m    i  rPw  dw 

J|  r-s  | 

=  ^  \  dv  sin  vr  sin  vs  \  dp  cos  vp  I  e~^   dw 
"0  Jo         J  p 

=  ^  \  dv  sin  vr  sin  vs  \  dw  sin  vw  . 


The  monotonicity  of  the  exponential  makes  the  inner  integral 
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positive  except  for   p  =  0  and   av  =  kit  ,  when  it  is  zero. 

Then 

„po         j  a.  .2    .2a 

(f,f*)    =  i|C2    l   ^   A     ds   0(s)    sin   vs  )     \    dw  e"Pw   sin  vw 

Jo       [yo  /     ^0 

>  0  for     p  >   0   . 

If     p   <   0   ,    then  let      p   =  -y      and   tne   kernel   of    the 
quadratic   form  is 

1    (eY|r+s|    .   eY|r-s|j    =  J£  slnh  yr  sinh  Ys   +  i   (e^l  r~s  I  -e^l  r+s  I  ) 

In  this    case 

(f,f*)    =^      Us^ls)    sinh  ys  ]     + 


oo         ,b.  2      2a 

i     ^/T'ds   0(s)    sin   vsj      V 


+c2    \     ^/(  ds   0(s)    sin   vs  j      \      dw  e"YW 


for     y  >   °   • 
Thus,    for  each      p.      the   characteristic   functions   of 
the   operator  and   its   adjoint   can  be    taken  to  be 
biorthogonal    ([10],    p„266). 

Theorem  3.8.      The    index  of    the    real   eigenvalues    of     A 
is    one. 

Proof.      If   the   index  of      p      Is    greater   than  one,    then  there 
is   an     f     such  that 

(p    -   A)2    f  =   0      and      (p    -   A)    f   j£  0 
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([11],  Th.  10.8),   Let   (3  have  multiplicity  s   and 


(p  -  A)  f.  =  0  for  j  =  1,2,0  0.,   s  >  1  ; 

then  also 

(P  -  A*)  f*  =  0   for  j  =  1,2,  ...,  s  >  1 

where   f"(x,  y)   is  a  linear  combination  of  the 

f .  (-x,y)  o   Since   (p  -  A)  f   is  a  characteristic  function 

of  A  , 


If   the   sets      [f .]      and      [f. ]      are   biorthogonal,    then 

J  J 


•&  .  AS- 


&.   =  ((p-A)f,fj)  =  (f,  (p-A")f\,)  =0   j  =  1,2, 
and   (p  -  A)  f  =  0  contrary  to  hypothesis» 
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l\.0      Remarks 

At  this  writing  it  has  been  impossible  to  prove  the 
existence  or  non-existence  of  complex  characteristic 
values  for  Re(\)  <  0  0   That  no  naive  extension  of  the 
method  used  by  Lehner  and  Wing  can  be  used  is  shown  as 
follows o   Let 


-N, 


h  =  \  +  "-x 


where      St       and     N.       are    real   symmetric   operators.      If 
for  real,    positive      p     and  for   real  functions      0     and     \J/ 


1^(0+140    =   p(j^i^) 


then 


pfi  =  R  0   -   N^\if 


and 


pt   =   N%0  +  M^    . 
Taking   inner  products   and   subtracting  yields 
(0,\fl)    +    (V,\V)    =   Im((j2f  +   iv|/)    ,    1^(0   +   ir|0)    =    0 


Therefore,    if     JL       is   definite    (as    it   is    for     Re(X)    >   0   , 
Im(\)    ^   0)    ,      there   can  be   no   real,    positive   characteristic 


3k 


value   p  o   This  is  the  essence  of  the  proof  given  by 
Lehner  and  Wing.   In  general,  however,   (0+i\J<,  L  (0+i\|f ) ) 
is 

(i)   an  entire  function  of  \ 

(ii)   not  a  polynomial 

(iii)   real  and  monotone  decreasing  for  \      real. 
The  function  must  assume  every  real,  positive  value  (with 
the  possible  exception  of  a  single  value)  infinitely  often. 
Since  it  can  assume  a  real  value   p   at  most  once  on  the 
real  line  and  not  at  all  for  the  right  half  plane  off  the 
real  line,  it  must  take  this  value  in  the  left  half  plane. 

It  is  possible  to  exclude  from  the  spectrum  a  band 
about  the  negative  real  line  (and  excluding  that)  whose 
width  decreases  exponentially  as   Re(X)  — >  -co.   By 
use  of  comparison  kernels  it  can  be  shown  that  the 
imaginary  part   It   of  the  transformation  is  definite 
in  this  band. 
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